Heat conductivity of the spin-Peierls compounds TiOCl and TiOBr 
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We report experimental results on the heat conductivity k of the S — 1/2 spin chain compounds 
TiOBr and TiOCl for temperatures 5K < T < 300 K and magnetic fields up to 14 T. Surprisingly, 
we find no evidence of a significant magnetic contribution to k, which is in stark contrast to recent 
results on S = 1/2 spin chain cuprates. Despite this unexpected result, the thus predominantly 
phononic heat conductivity of these spin-Peierls compounds exhibits a very unusual behavior. In 
particular, we observe strong anomalies at the phase transitions T c i and T C 2- Moreover, we find an 
overall but anisotropic suppression of k in the intermediate phase which extends even to temperatures 
higher than T C 2. An external magnetic field causes a slight downshift of the transition at T c \ and 
enhances the suppression of « up to T C 2- We interprete our findings in terms of strong spin-phonon 
coupling and phonon scattering arising from spin-driven lattice distortions. 

PACS numbers: 66.70.-f, 75.40.Gb, 75.10.Pq, 68.65.-k 
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I. INTRODUCTION 



Understanding low-dimensional quantum spin- 1/2 sys- 
tems is one of the challenges of contemporary condensed 
matter physics. In particular, transition metal oxides 
provide a rich playground for studying novel phenom- 
ena, arising from the interplay between lattice, orbital, 
spin, and charge degrees of freedom. The recent dis- 
covery of a substantial magnetic heat conductivity K mag 
in ID quantum spin systems 1 - - — together with the the- 
oretical prediction of ballistic transport in ID S = 1/2 
Heisenberg chains 1 ^— has caused intense experimental 
and theoretical research on the behavior of these sys- 
tems. The best experimental realizations of S = 1/2 
systems showing magnetic heat transport are up to now 
found among copper-oxides (cuprates) such as the spin 
chains SrCu02 and S^CuOa 2 ' 11 and the spin ladder com- 
pounds (Ca,La,Sr)i4Cu2404if 1 -'££'£ Characteristic for the 
cuprate chain systems is a Cu3c£ 9 configuration which 
gives rise to S = 1/2 and a large exchange coupling 
J/fc B w 2000 K along the chain/ladder direction. As 
a consequence of this quasi ID magnetic structure, these 
systems exhibit a strongly anisotropic thermal transport 
behavior. Perpendicular to the low-dimensional spin 
structure a typical phononic heat conductivity K p h is 
found. However, parallel to the low-dimensional spin 
structure the heat conductivity is strongly enhanced even 
up to room temperature since a large K mag adds to K p h- 
These in many aspects excellent realizations of S = 1/2 
Heisenberg chains do not undergo a spin-Peierls transi- 
tion, i.e. a transition to a spin-dimerized ground state at 
the expense of a lattice distortion that normally should 
arise from the spin-phonon coupling of a spin chain and 
the phonon system in which it is embedded. Surprisingly, 
only one Cu based spin system, CuGeOa, is known to ex- 
hibit a spin-Peierls transition^. The exchange energy of 
this compound is J/fcs ~ 160 K 16 and the transition 



to the non- magnetic ground state is at T c « 14 K. The 
heat conductivity of CuGeOa has been studied by several 
groups 1 ^— with controversial results. ID magnetic heat 
conductivity has been suggested to give rise to a signif- 
icantly enhanced heat conductivity at T < T c 17 . How- 
ever, the observed low-temperature peak has been shown 
to be present both in the heat conductivity parallel and 
perpendicular to the chain and thus can be rationalized 
in terms of phononic transport alone ^ 

Also spin 5 = 1/2 systems, but based on early transi- 
tion metal ions with electronic configuration 3d 1 , the tita- 
nium oxyhalides TiOX, with X=Br or CI shifted recently 
into focus. These compounds are considered as good re- 
alizations of S — 1/2 spin chains which are formed by di- 
rect overlap of Ti t2 g orbitals along the crystallographic 
b direction 2 ^— with rather high magnetic exchange cou- 
pling J (CI) w 676 K22i2i an d j(B r ) w 375 K 23 ~ 25 . The 
compounds undergo two phase transitions T c i, T C 2 23 i 26 
where the lower one, at T c i, leads to a non-magnetic 
dimerized state 20 which is accompanied by a doubling of 
the unit celL 24 ' 27 These features thus render the Ti oxy- 
halides the second (besides CuGeOa) type of inorganic 
compounds which undergoes a spin-Peierls transition. 
However, as compared to CuGeOa the dimerized state 
occurs at much higher temperatures, viz. T c i a = 67 K 
for TiOCl and T cljBr = 28 K for TiOBr. However, sev- 
eral experimental results are inconsistent with a canoni- 
cal spin-Peierls scenario. There are two successive phase 
transitions and the transition to the non-magnetic state 
at T c i is of firs t 24 ' 26 ! 27 and not of second order as in 
CuGeOa. Interestingly, in the intermediate regime be- 
tween T c i and T c i an incommensurate superstructure is 
found. 28 - 31 Above T c2 (T c2:C i = 91 K for TiOCl and 
7c2,Br = 48 K for TiOBr) the system is in a pseudo 
spin-gap regime up to a characteristic temperature T* 
which for TiOCl extends up to T* » 135 K with a 
large singlet-triplet energy gap of E g = 430 Ki 30 ' 32 " — 
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Figure 1: (color online). Thermal conductivities n a (o) and Kb (♦) of TiOCl and TiOBr as a function of T. The insets depict the 
behavior of the thermal conductivity around the respective phase transitions. Insets (a) and (c) show the hysteretic behavior 
around the phase transition at T c \ which is characteristic for a first-order transition. The arrows mark the corresponding curves 
for cooling and heating. The lower insets (b) and (d) show k, around the transition at T C 2 also measured for cooling and heating. 
It is without hysteresis and therefore the transition is of second order. 



First explanations of the intermediate phase proposed 
orbital fluctuations but this has been ruled out by op- 
tical measurements in combination with cluster calcula- 
tions that showed, that the crystal field splitting is large 
enough to quench the orbital degree of freedom. 23 Re- 
cent explanations focus on the interplay between intra- 
and interchain frustrations and a related dimensionality 



II. EXPERIMENT 



crossover 
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The relatively high magnetic exchange constants of the 
titanium oxyhalides render them good non-cuprate can- 
didates for exhibiting a sizeable magnetic heat conduc- 
tivity arising from the ID S = 1/2 spin chains. In this 
paper, we experimentally investigate the thermal con- 
ductivity k of TiOCl and TiOBr with a special focus 
on potentially arising magnetic contributions to k. Sur- 
prisingly, no indication for magnetic heat transport is 
observed and we find instead that k is dominated by 
phononic heat conduction. However, strong anomalies 
occur at the phase transitions T cl and T C 2 and we find an 
overall suppression of the phononic k which is anisotropic 
in the incommensurate phase and which extends to tem- 
peratures higher than T c i. For TiOBr the application of 
an external magnetic field of 14 T slightly shifts T c \ to- 
wards lower temperature and causes a weak further sup- 
pression of k in the intermediate regime. 



Single crystals of TiOCl and TiOBr were synthesized 
by a chemical vapor transport technique leading to small 
plate-like crystals^ The crystallinity was checked by x- 
ray diffraction. Typical crystal dimensions are a few 
mm 2 in the afe-plane but only around 20 /im along the 
c axis. Rectangular samples with typical dimensions of 
(2 x 1 x 0.02) mm 3 with the longest side being parallel 
to the a and b axis, respectively, were cut from the crys- 
tal plates. Measurements of the thermal conductivity 
as a function of temperature T in the range of 7-300 K 
were performed with a standard four probe technique^. 
Because of the small thickness of the crystals the usual 
uncertainty of 10% for k due to the error in the determi- 
nation of the crystal geometry is exceeded by some ex- 
tent. Furthermore, the small thickness along the c axis 
also prevented to measure k along this direction. In or- 
der to compare the anisotropy of k along the a and b 
directions the individual samples were cut from the same 
crystal plate thus keeping the relative error between the 
two directions small. The mounting of TiOBr into the 
heat conductivity probe was performed under Argon at- 
mosphere in order to minimize degradation of the sample. 



3 



III. RESULTS 

Figure [T] shows the temperature dependence of the 
thermal conductivities along the a and b and Kb) 

of TiOCl and TiOBr in zero magnetic field. We focus first 
on the results for TiOCl which are shown in the left panel 
of Fig.[TJ A first glance at the data already suggests that 
the temperature dependence of n is governed by the two 
phase transitions at T c \ and T C 2 which divide the data 
into three regimes. At low temperature the heat conduc- 
tivity parallel to the chains, Kb, exhibits a strong peak 
at ~ 25 K with a maximum value Kb ~ 58Wm _1 K _1 
which is a typical feature of a phonemic heat conductiv- 
ity Kph at low temperature. It arises from two competing 
effects: At very low temperature the mean free path of 
phonons is determined by the crystal boundaries and de- 
fects and therefore is practically T- independent. Hence, 
K p h increases due to the increasing number of phonons. 
At higher temperature the mean free path is T-dependent 
as the number of umklapp processes rises exponentially. 
This overcompensates the effect of a rising phonon pop- 
ulation and thus K p h decreases with further rising T, i.e., 
K p h usually shows a maximum. Interestingly, Kb deviates 
from this conventional behavior at T c i, where a sharp 
drop occurs to about 60% of the value of Kb at just below 
the transition. In the intermediate phase Kb continuously 
decreases further with rising T just until T C 2 is reached. 
Upon rising T through T c i we find that Kb changes slope 
and exhibits a weak increase in the entire high tempera- 
ture phase, i.e., at T > T C 2, up to room temperature. 

A very similar temperature dependence is observed in 
the heat conductivity perpendicular to the chains, K a . In 
this case, the peak at T w 25 K is somewhat smaller 
(«o,max ~ 48Wm _1 K -1 ) than that in Kb- A similarly 
sharp drop as in the latter occurs at T cl . However, the 
actual drop at the transition is relatively weaker as in 
the other direction. Interestingly, despite a similar fur- 
ther decrease of K a when rising T towards T C 2 as in Kb, 
we find that K a remains always somewhat larger in this 
intermediate regime. The slope of K a changes at T C 2, but 
remains negative up to T rs 150 K, in contrast to the 
findings for Kb (cf. Fig. [TJb)). 

Before discussing these pecularities in detail, we briefly 
summarize the results for TiOBr which are shown in the 
right panel of Fig.[T] The general T-dependence of k has 
large similarities with that of TiOCl, including the ob- 
served anomalies. There are, however, slight differences 
which are worth to be pointed out: First, the phononic 
peak of both K a and Kb of TiOBr is by a factor of about 
4 larger than that in TiOCl and is located at some- 
what lower temperature (~ 17 K). Both features point 
to a lower defect density in the case of TiOBr. This 
is corroborated by room temperature x-ray diffraction 
which showed much sharper spots for TiOBr. Second, 
at T < T C 2 the anisotropy between K a and Kb is similar 
to that of TiOCl. More specifically, at T < T cl we find 
K a < Kb, and K a > Kb at T c \ < T < T C 2, i.e. the drop at 
T c i and the reduction of k are relatively stronger in Kb 



than that in K a . Interestingly, the anomaly in Kb at T C 2 
is much stronger than that in TiOCl since a clear dip is 
observable at the transition (cf. Fig.QJd)). Moreover, in 
contrast to TiOCl we observe that both K a and Kb de- 
crease with rising temperature at T > T C 2 up to room 
temperature where K a remains slightly larger than Kb- 



IV. DISCUSSION 

The overall very weak anisotropy of the k data sug- 
gests without further analysis the unexpected conclusion 
that magnetic heat transport in the spin chains of this 
material is negligible in both TiOCl and TiOBr. Oth- 
erwise a significant enhacement of Kb with respect to K a 
should occur since heat transport by magnetic excitations 
is only expected along the ID spin chain, i.e. parallel to 
b. One might speculate that the weak anisotropy that 
is present in the low temperature regime T < T c \ is the 
indication of a weak magnetic contribution along b which 
could give rise to the observed Kb > K a . However, the ob- 
served anisotropy by a factor ~ 1.2 matches that of other 
phononic heat conductor s 40 ' 41 and can conventionally be 
explained by differences in the phonon velocity. 

At higher temperatures (T > T c i) magnetic contri- 
butions appear even more unlikely, since in all cases 
K b iS Ka- However, in this regime a small magnetic 
contribution to Kb might still be present if the expected 
anisotropy was masked by differences in the phononic 
transport along the two crystallographic directions. Con- 
centrating only on the thermal conductivity Kb we esti- 
mate the thus maximum possible K mag by performing a 
phononic fit based on the so-called Callaway model 4 ^ to 
the low temperature peak and extrapolate this fit towards 
room temperature. The fit is depicted by the solid line 
in Fig. [1] and yields a very good agreement up to T c i but 
deviates strongly from the data at higher temperatures. 
In particular, at high temperatures (T > 180 K) the fit 
is clearly lower than the data. We use the difference be- 
tween the fit K p h,Fit and the data at room temperature to 
obtain an upper estimate for the possible magnetic con- 
tributions K mag = Kb — K p h,Fit- In order to analyse the 
thermal transport we estimate the magnetic mean free 
path Z mag using an approximation of K mag of a S = 1/2 
Heisenberg chain^ii - 

_ 2n s fc B 2 ; f^ST 2 exp(a;) 

mas " nh mas Jo X (exp(x) + l) 2 ' l j 

where n s is a geometrical factor that counts the num- 
ber of chains per unit area. For both compounds 
this yields a negligibly small mean free path of only 
2-3 lattice constants. 53 Considering the fact that the 
Callaway model usually underestimates k p ^ at room 
temperatur e) 43 ' 44 and that Kb < K a at higher temperature 
any realistic value for the mean free path should be even 
smaller which essentially rules out magnetic transport in 
the Ti oxyhalides. 
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Figure 2: (color online). Thermal conductivity func- 
tion of T in TiOBr with (nb.uT (A)) and without («;{, (♦)) an 
applied magnetic field of 14T along the chain direction. Inset 
(a) illustrates the shift of T c i towards lower temperatures in 
the presence of a magnetic field. In inset (b) A/t = kj, — Kb,iiT 
shows the decreasing influence of the magnetic field on Kb in 
the intermediate regime. The curves used in the substraction 
are from the measurements that approach the phase transi- 
tions from low temperatures. 



There are not many scenarios which straightforwardly 
explain this unexpected result. The absence of magnetic 
heat conduction in magnetic materials has been discussed 
by Sanders and Walton in terms of a very large magnon- 
phonon relaxation time 45 . It is obvious that this sit- 
uation cannot be realized in Ti-oxyhalides since a sig- 
nificant spin-phonon coupling must be present in these 
compounds to allow for a spin-Peierls transition at con- 
siderably high temperatures. In fact, it is therefore more 
reasonable to explain the absence of magnetic heat con- 
duction by a particularly strong spin-phonon coupling 
which gives rise to strong scattering of spin excitations 
and thus prevents the magnetic heat conduction. One 
might speculate that even more exotic excitations such 
as orbital fluctuations are relevant for suppressing K mag . 
We point out, however, that orbital excitations have been 
shown to be unimportant for the low-energy physics in 
these compounds! 22 ' 23 

The negligible magnetic heat conduction in the Ti- 
oxyhalides implies that the unusual temperature depen- 
dence and also the slight anisotropy should be rational- 
ized in terms of pure phonon heat conduction, which 
has been proven to be a sensitive probe to peculari- 
ties of the lattice such as superstructures and disor- 
der i 40 i 41 ' 46 i 47 The considerable jump in k at T c \ clearly 
indicates that the phonon heat conduction in the in- 
termediate phase is strongly suppressed with respect to 
that of the commensurate dimerized phase at T < T c \ 



where ordinary phonon heat conduction is observed. 
This reflects the abrupt transition towards a lattice with 
strongly disturbed periodicity and anharmonicity which 
causes enhanced phonon scattering and is entirely con- 
sistent with the incommensurate lattice distortion in this 
regime! 29 ' 32 ' 48 We have investigated the nature of this 
phase transition at T c \ and find for both compounds a 
clear hysteretic behavior which confirms the transition 
being of first order (see Fig. [TJa) and[IJc)). Such first- 
order character has already been reported from magnetic 
susceptibility, specific heat, thermal expansion and x-ray 
data of the superstructure satellites! 22 ' 23 ' 26 ' 27 Since the 
magnetic exchange is smallest in TiOBr we have searched 
for possible effects of a magnetic field on K b . As is de- 
picted in Fig. [2] a magnetic field of B = 14 T along the 
b direction has only little influence on the thermal con- 
ductivity K6 ; i4T- However, we detect a slight downshift 
of the phase transition at T c \ by ~ 400 mK which is con- 
sistent with a downshift of ~ 130 mK that has been re- 
ported from x-ray diffraction at B = 10 T for TiOCl. 49 
Moreover, starting at T c i, Kb.iiT is slightly smaller com- 
pared to the measurement without field, but gradually 
approaches it for increasing temperature. In Fig. &h) 
the difference Ak = K b — Kb,iiT between both curves is 
shown, illustrating the decreasing influence of the mag- 
netic field, until it vanishes at T C 2- This suggests that the 
spin-induced incommensurate lattice distortion in this in- 
termediate phase is further enhanced by an external mag- 
netic field. 

The thermal conductivity across the phase transitions 
at T C 2 shown in more detail in Fig. [Ijb) andQTd), does 
not exhibit a hysteretic behavior which is indicative of a 
second-order transition. The overall impact of this tran- 
sition on k is much smaller than that at T c \. Interest- 
ingly, in the high-temperature phase above T c i the ther- 
mal conductivity appears still significantly suppressed 
with respect to the low-temperature phase at T < T c ±. 
In Fig. [T] this is clearly seen when comparing the data to 
the phononic fit which remains much larger than k up 
to T* ~ 100 K and T* ~ 180 K for TiOBr and TiOCl, 
respectively. Only at higher temperatures a more typi- 
cal behavior is observed with Kph.Fit < k< The apparent 
suppression of k in the regime T C 2 < T < T* clearly in- 
dicates, that strong phonon scattering occurs despite the 
absence of any static long range lattice distortions. A 
reasonable origin of this enhanced scattering are precur- 
sors of the spin-Peierls transition, either as short-range 
static lattice distortions or as slowly fluctuating precur- 
sors (soft phonon type) . This is consistent with the pseu- 
dogap seen in magnetic resonance measurement a 21 ' 32 ' 33 
and incommensurate structural fluctuations found by x- 
ray diffraction^ 8 - 

In both compounds the suppression of k in the inter- 
mediate phase is clearly anisotropic, since the drop of Kb 
at T c i is relatively stronger as compared to n a and Kb<K a 
in the entire phase where Kb of TiOBr even shows a local 
minimum at T C 2- Similar anisotropic scattering has pre- 
viously been observed, e.g., in stripe-ordering compounds 
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which possess anisotropic correlation lengths of the stripe 
order close to the transition^ The stronger suppression 
of Kb than K a in the present case can be understood by 
looking at the modulation amplitudes for TiOCl^ and 
TiOBr in the incommensurate phased. Those indicate 
that the shifts of the atoms out of the periodic position 
of the structure at room temperature are larger in the 
direction of the b axis than those along the a axis. The 
resulting larger anharmonicity along b is likely causing in- 
creased scattering and therefore the observed lower ther- 
mal conductivity. 

There is a slight difference in the thermal conductivity 
between both compounds near room temperature where 
phonon scattering arising from the spin-Peierls transi- 
tion can be considered to be relatively weak. For TiOBr 
the slope of k is negative while it is positive for TiOCl. 
At the same time the absolute value of k is significantly 
higher in TiOBr. This corroborates the previous conclu- 
sion that our TiOBr crystals have a lower defect density 
than the TiOCl ones because the observed temperature 
dependence of k for TiOBr is much closer to the expected 
oc T^ 1 decrease of a clean phonon heat conductor^ On 
the other hand, the lower k of TiOCl with a weak posi- 
tive slope is typical for more disordered heat conductors, 
where also rather small contributions to k, such as heat 
transport by optical phonons^i become relevant. 

V. SUMMARY 

In conclusion, we have shown that the magnetic ther- 
mal conductivity in the TiOX is negligible due to strong 



spin-phonon scattering. The heat transport can thus 
be understood in terms of pure phononic conductivity. 
At the phase transitions we find strong anomalies which 
are consistent with the lattice distortions. Starting at 
low temperatures, the first phase transition T c % towards 
the dimerized state can be shifted to lower tempera- 
tures by an external magnetic field. Additionally, this 
leads to a slight suppression of the thermal conductivity 
throughout the intermediate regime and gradually gets 
smaller when approaching T c2 . Comparing the measure- 
ments along the different crystallographic directions in 
this regime, the stronger suppression along Kb for both 
compounds is consistent with a higher incommensurabil- 
ity of the lattice in this direction. Finally, by a com- 
parison of the extrapolated thermal conductivity from a 
phononic model to the measurement at higher tempera- 
tures it was argued that the thermal conductivity is still 
supressed up to a temperature T* which is either a sign 
of short-range lattice distortions or phonon softening. 
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